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Abstract 

The fundamental assumption of conventional transition state theory is the 
existence of a dividing surface having the property that trajectories origi- 
nating in reactants must cross the surface only once and then proceed to 
products. Recently it has been shown ^12] how to construct a dividing sur- 
face in phase space for Hamiltonian systems with an arbitrary (but finite) 
number of degrees of freedom having the property that trajectories only 
cross once locally. In this letter we provide an argument showing that the 
flux across this dividing surface is a minimum with respect to certain types 
of variations of the dividing surface. In this sense, conventional transition 
state theory is shown to be equivalent to variational transition state theory. 

PACS: 82.20.Db, 82.20.Nk, 05.45.- 



1 Introduction 

Transition state theory provides a fundamental framework for computing 
chemical reaction rates. The original ideas are due to Wigner, Polanyi, 
and Eyring, yet much work on different aspects of the subject continue to 
this day (see the recent reviews of [3] and |H). In recent years transition 
state theory has been shown to be much more widely applicable than just 
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for problems in chemical reactions. It has been used in atomic physics [3], 
studies of the rearrangements of clusters jS] , solid state and semi-conductor 
physics 0|H|) cosmology 9 , and celestial mechanics |1U| . 

The fundamental assumption of transition state theory is the existence 
of a dividing surface having the property that trajectories originating in 
reactants must cross the surface only once and then proceed to products. 
The reaction rate is proportional to the flux through the dividing surface. 
The construction of such a surface in specific systems, especially those with 
more than two degrees of freedom (DOF), has posed extreme difficulties. 
The phenomenon of recrossing of a dividing surface gives rise to a larger value 
for the flux, which naturally inspires the idea of varying the dividing surface 
so that the flux is minimised, this has lead to the subject of variational 
transition state theory ( JTJ, see also the review of [T2]). 

For systems with 2 DOF described by a Hamiltonian H of the simple type 
kinetic-plus-potential the problem of constructing a dividing surface (transi- 
tion state) with trajectories only crossing once and having minimum flux was 
solved during the 70's by McLafferty, Pechukas and Pollak [T31 ITU IT51 Hfl] . 
They considered the line in configuration space given by the projection of 
an unstable periodic orbit (the Lyapunov periodic orbit associated with a 
saddle point), see figure UJ^a). This line which is referred to as periodic orbit 
dividing surface (PODS) separates the reactant region (x < 0) from the 
product region (x > 0) by connecting two pieces of an equipotential. This 
solves the problem of (local) recrossings as there is no trajectory evolving 
from reactants to products (or vice versa) whose projection to configuration 
space is tangent to this line. After a reacting trajectory has crossed the 
PODS it has to leave the neighbourhood of the PODS before it can possibly 
recross it. 

To comprehend a PODS as a phase space object we rewrite the energy 
equation as 

p 2 x + P l = 2m(E-V(x,y)). (1) 

For (x, y) on the projection of the periodic orbit in figure^a) the equation in 
(fT| defines a family of circles in (p x ,p y ) which shrink to points at the turning 
points of the periodic orbit, i.e. the PODS is a 2-sphere. The periodic orbit 
which has p x = can be considered as the equator of the PODS. It divides 
the PODS into two hemispheres which have p x > and p x < 0, respectively. 
The flux from x < ("reactants") to x > ("products") through the PODS 
in figure^ a) is commonly written as 

T E = I dxdydp x d Py 5{E - H)5(x - x o )0fe)— • (2) 



2 




Figure 1: (a) Equipotentials and PODS (bold line) near a saddle, (b) Saddle 
plane (<?i,pi) with projections of the (2n — 1) -dimensional energy surface 
H = E for E > (all but the hatched region), the stable and unstable 
spherical cylinder (the pi-axis and the gi-axis), the NHIM (the origin), and 
the (2n — 2)-spheres p\ — q\ = (dividing surface). The light and dark 
grey regions are the projections of the energy surface volumes enclosed by 
the forward and backward reactive spherical cylinder Wf(E) and Wb(E), 
respectively, pi—qi = ±c, c > 0, delimit the region of validity of the normal 
form. 



The step function @(p x ) effectively restricts the integral to one hemisphere 
of the PODS. Using Stokes' theorem it is not difficult to see that the flux in 
PJ) is given by the action of the periodic orbit (p.o.), i.e. 




pdq. 



Unfortunately, the methods of McLafferty, Pechukas and Pollak are spe- 
cific to 2 DOF, and as well to Hamiltonian's of a certain form. Recently 
it has been shown ^ |2| how to construct a dividing surface in phase space 
for Hamiltonian systems with an arbitrary (but finite) number of DOF hav- 
ing the property that trajectories only cross once locally. A computational 
algorithm for realising this construction was also given. However, it was 
not shown that this transition state was a surface of minimum flux. In this 
letter we give a geometrical argument demonstrating that this is indeed the 
case, and in the process show that for the phase space transition states con- 
structed according to the techniques in ^ [2| , conventional transition state 
theory is equivalent to variational transition state theory. In order to do 
this we must first recall the elements of the many DOF theory in ^HJ. 
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2 Dividing surface for multi-dimensional systems 



In the general case of n DOF we consider an equilibrium point of saddle- 
center-... -center type, i.e. the linearised vector field about the equilibrium 
point has one pair of real eigenvalues ±A and n — 1 imaginary eigenval- 
ues ±iuJk, k = 2, . . . ,n, where without restriction of generality A, u>k > 0. 
Equilibria of this type are characteristic for systems with a reaction type 
dynamics, and they occur in all of the examples given in the introduction. 
The transport is controlled by various high-dimensional manifolds which can 
be realised and computed through a procedure based on Poincare-Birkhoff 
normalisation. In fact, under general conditions (see [2] for the details), near 
a saddle-center-... -center one can construct a sequence of local, non-linear, 
symplectic transformations of the phase space coordinates that transform 
the Hamiltonian into 

n 

H = \p 1 q 1 + 'Y J —(pl + ql) + f 1 (l,q 2 ,...,p n ,q2,---,qn) 

k=2 

+ f2(q2,---,p n ,q2,---,q n ) (3) 

up to any desired order. Here (gi, . . . , q n ,Pi, ■ ■ ■ ,Pn) are canonical phase 
space coordinates, X = p\q\ and /i, f 2 are at least of third order, i.e. 
they are responsible for the nonlinear term in the Hamiltonian vector field. 
Moreover f\ has the property that it vanishes for 2 = 0. (qi,p\) play the 
role of reaction coordinates. (q 2 , ■ ■ ■ ,q n ,P2, ■ ■ ■ ,Pn) are the bath coordinates. 

The importance of saddle-center-... -center equilibria for reaction type 
dynamics can be inferred from the topology of energy surfaces H = E. 
First consider the quadratic Hamiltonian given by the first part in (jSJ) and 
write the energy equation as 

E + j(Pi ~ H) 2 = \(Pi + H) 2 + E f (pl + Qi) ■ (4) 

k=2 

For E < the left hand side is positive for p% — q% < — (— AE/X) 1 / 2 or 
Pi ~ qi > ( — 4.E/A) 1 / 2 . For a fixed p\ — q\ in either of these ranges the 
equation in (JIJ defines a (2n — 2) dimensional sphere S 2n ~ 2 . For p\ — q\ = 
db(— AE/X) 1 / 2 the (2n — 2)-spheres shrink to points. The energy surface thus 
appears as two disjoint spherical cones which correspond to "reactants" and 
"products", respectively. Increasing the energy to E > 0, the left hand 
side of Q is strictly positive. The formerly disjoint components merge 
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and the energy surface becomes a spherical cylinder S 2n ~ 2 x R. Restricting 
to a sufficiently small neighbourhood by confining p\ — q± to an interval 
/ = [— c, c] with c > sufficiently small and with E sufficiently close to zero, 
the topological consideration remain true if the non-linear terms f\ and fi 
are taken into account. Moreover, for a high but finite order normal form, 
the error from neglecting the non-normalised "tail" of the Hamiltonian can 
be made as small as desired by choosing the interval [— c, c] sufficiently small 
and E sufficiently close to zero. 

The importance of the normal form arises from the fact that it gives 
explicit expressions for all the manifolds which control the dynamics from 
reactants to products. For a fixed energy E > these manifolds are 

• The saddle sphere Sf^fj^E): On the energy surface H = E the 
equation p\ = q\ = defines a (2n — 3)-sphere which we denote by 
^nhimC^)- I* can ^ e considered as a "big saddle". In fact, it is a so- 
called normally hyperbolic invariant manifold (NHIM) where normal 
hyperbolicity means that the expansion and contraction rates trans- 
verse to the manifold dominate those tangent to the manifold. 

• The forward and backward reactive spherical cylinders Wf(E) and 
Wb(E): The saddle sphere has stable and unstable manifolds W S (E) 
and W U (E) which are iso-energetic, i.e. contained in the energy sur- 
face, and which are explicitly given by q\ = and p\ = 0, respec- 
tively. These invariant manifolds have the topology of spherical cylin- 
ders S" 2 " --3 x /. Since they are of codimension 1 in the energy sur- 
face, i.e. they are of one dimension less than the energy surface, they 
act as impenetrable barriers. W S (E) and W U {E) each appear as two 
branches. We call the branch of W S {E) which has p\ > the forward 
branch Wj(E) and the branch which has p\ < the backward branch 
W b s {E). Likewise the forward branch Wj(E) of W U (E) has qi > 
and the backward branch W b u {E) has q\ < 0. We call the union of 
the forward branches, Wf(E) := WJ(E) U Wj(E), the forward reac- 
tive spherical cylinder. Similarly, we call the union of the backward 
branches, Wb(E) := W b s {E) U W b u (E), the backward reactive spherical 
cylinder. The significance of these spherical cylinders arises from the 
fact that they enclose volumes of the energy surface which contain 
all forward and all backward reactive trajectories, respectively. All 
non-reactive trajectories are contained in the complement of these two 
volumes. 

We define the dividing surface or transition state as follows. 
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• Transition state S 2 ™~ 2 (E): On the energy surface H = E the equation 
pi = q\ defines a (2n — 2)-sphere which we denote by S 2 ™~ 2 (E). It is 
of codimension 1 in the energy surface. It divides the energy surface 
into two components: the reactant region p\ — q\ > and the product 
region pi — q\ < 0. The saddle sphere ^nhim^) can ^ e considered 
as the "equator" of the transition state. It divides S 2 ™~ 2 (E) into 
the two hemispheres: the forward hemisphere D 2 ™J 2 (E), which has 
Pi = qi > 0, and the backward hemisphere D 2 ™^ 2 (E), which has p\ = 
qi < 0. D^J 2 (E) and D 2 ^ 2 (E) are topological (2ra-2)-discs. Except 
for its equator, which is an invariant manifold, the transition state is 
everywhere transverse to the Hamiltonian flow as is easily seen from 
the equations of motions derived from the normal form Hamiltonian 
©. This means that a trajectory, after having crossed the transition 
state, has to leave the neighbourhood of the transition state before it 
can possibly cross it again, i.e. the transition state locally is a "surface 
of no return" . 

An important advantage of the normal form coordinates is that dynam- 
ical issues related to flux and reactivity can be understood to a great extent 
from the projection to the plane of the reaction coordinates (qi,pi), see 
figure Wlb). Due to the Constance of the saddle integral 2 = piqi trajec- 
tories project to hyperbolas. Forward and backward reactive trajectories 
project to the first quadrant qi,Pi > 0. It approaches the forward hemi- 
sphere D 2 ™J 2 (E) of the transition state inside of the forward branch of the 
stable cylinder Wj(E) whose interior projects to the part of the first quad- 
rant above the diagonal p\ = q\ > 0. After crossing the transition state 
the trajectory leaves inside of the forward branch of the unstable spherical 
cylinder WJ\E) whose interior projects to the part of the first quadrant 
below the diagonal p\ = q\ > 0. While similar considerations hold for 
the backward reactive trajectories, non-reactive trajectories project to the 
second quadrant p\ > 0, q\ < which corresponds to reactants or to the 
fourth quadrant p% < 0, q\ > which corresponds to products. Therefore 
reactive trajectories have T = p\q\ > and non-reactive trajectories have 
I = PiQi < 0. 

3 Minimal flux property of the transition state 

Consider at first the phase space volume form Q = dpi A dqi A • • • A dp n A dq n , 
which in terms of the symplectic 2-from uj = Yl"k=i dPfc^d^fc can be written as 
= uj n /n\. Let j] be an energy surface volume form defined via the property 
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dH A 77 = $1. Then the flux through a codimension 1 submanifold of the 
(2n — l)-dimensional energy surface H = E is obtained from integrating 
over it the "flux" form £1' given by the interior product of the Hamiltonian 
vector field Xh with 77 ^7j) i- e - 

n' = i XHV = j^jy^- 1 (5) 

where ix H r](£i, . . . , £2^-2) = 7/(fi, . . . ,^2u-2,Xh) for any 2n - 2 vectors 
The second equality in (J5J) is easily established on a non-critical energy 
surface, i.e. on an energy surface which contains no equilibria. The flux 
form Q' is exact. In fact the generalised "action" form 

n 1 
<P = y^Pkdqk A TT7 W "~ 2 

^-^ (n — 1 

fe=l v ' 

has the property dip = Q! and facilitates the use of Stokes' theorem to 
compute the flux. In the case of two DOF we simply have Q! = u) = dpi A 
dqi+dp2 Adq2 and <p = "P\dq\+ , p2dq2 • Since the transition state S 2 ™~ 2 (E) is a 
sphere, that is, a manifold without boundary, it follows from Stokes' theorem 
that the integral of O' over S 2 ™~ 2 (E) is zero. As in the case of PODS 
one has to distinguish between the directions in which the Hamiltonian 
flow crosses the transition state. Given a normal bundle over S 2 ™~ 2 (E) 
the direction can be specified by the sign of the scalar product between 
the normal vectors and the Hamiltonian vector field. This scalar product 
is strictly positive on one hemispheres of S 2 ™~ 2 (E), strictly negative on the 
other hemisphere and zero only at the equator of S 2 ™~ 2 (E), i.e. at the saddle 
sphere S 2 J^ A (E), where the Hamiltonian vector field is tangent to S 2 ™~ 2 (E). 
Likewise, the flux form f2' on S 2 ™~ 2 (E) vanishes nowhere on D 2 ™7 2 (E) and 

D tsb 2 ( E ) and is identically zero on S^him (-^) • ^ * s na t ura l to take as the 
orientation of D 2 ™7 2 (E) and D 2 ™7^ 2 (E) the orientation they inherit from the 
transition state. Without restriction we may assume that the orientation of 
S t ™~ (E) is such that Q' is positive on the forward hemisphere D t ™7 (E) 
and negative on the backward hemisphere D?™7 2 (E), i.e. Q 1 and — O' can 
be considered as volume forms on D 2 ™J 2 (E) and D 2 S \ 2 (£), respectively. It 
follows from Stokes' theorem that the flux through the forward and backward 
hemispheres, J D 2n-2 (E) Q' and f D 2n-2 (E - ) £1', have the same magnitude but 

ts, f >■ > ts.b > 

opposite sign and can be computed from integrating the action form <p over 

the saddle sphere: f n 2n-2,^ O,' = — f n 2n-2, m fl' = I f Q 2 n -3 fTr ^ <p\. We call 
JD ts,i ( h > JU ts.b W ^nhimW 
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the positive quantity J D 2n-2^ Q' the forward flux and the negative quantity 

J D 2n-2^ the backward flux through S 2 ™~ 2 (E) . 

We now show that the forward flux through S 2 ™~ 2 (E) is minimal. This 
can be stated as a variational problem and a result in this direction is ob- 
tained by MacKay JH] who proves that the integral J c ip of the action form 
<p over codimension 2 submanifolds C of the energy surface is stationary with 
respect to variations of C if and only if C is invariant under the Hamiltonian 
flow. Since the saddle sphere S'nhim^) ^ s an i nvar i an t manifold, MacKay 's 
result implies that the "action" of S^^ /L (E), f S 2ns ^ cp, is stationary with 

respect to variations of <5nhim(-^)- MacKay considers arbitrary variations 
within the energy surface and this leads to an indefinite variational princi- 
ple. In fact, in the context of transition state theory it is more useful to 
consider variations of the codimension 1 transition state S 2 ™~ 2 (E), which in 
a sense that will become clear below, imply variations of the codimension 2 
saddle sphere 5'nhim(-^') m a sma Uer, more suitable class than in MacKay's 
case. 

Consider at first the case of 2 DOF for which it is possible to visu- 
alise the energy surface S 2 x / in 3-dimensional space as a nested set of 
2-spheres parametrised along the interval / in radial direction. This is the 
so-called McGehee representation which is shown in the first panel of 
figure I2J We consider a slight generic, iso- energetic deformation of the tran- 
sition state S 2 S (E) that does not "preserve" the saddle sphere S 1 NHIM {E), 
i.e., Snhim^) is n °t entirely contained in the deformed transition state, see 
Figures 0c)- (e). The deformation, which we denote by S 2 S (E), is described 
mathematically as a graph over the transition state S 2 S (E), and therefore 
S 2 S (E) inherits the orientation of S 2 S (E) (which is important because we 
want to preserve certain aspects of the directionality of the flux, to be de- 
scribed below). Because S 2 S (E) is chosen to lie in the energy surface, it 
still separates the energy surface into two disjoint components like S 2 S (E). 
Also like S 2 S (E), the deformation S 2 S (E) contains a circle S l (E) of points at 
which the Hamiltonian vector field is tangent to S 2 S (E) . This can be proven 
analytically, but intuitively it is easy to see from figure |2^d). 5 ,1 (£') can 
be considered as the equator of S 2 S (E); it divides S 2 S (E) into forward and 
backward hemispheres D 2 S { (E) and D 2 S ^(E) on which the flux form Q' is 
strictly positive and negative, respectively. ^(E) can be considered as the 
deformation of 5'nhim(^) induced by the deformation of S 2 S (E). It should 
be realised that in contrast to the case of S^ mM (E), the Hamiltonian vector 
field is not tangent to the deformation S' 1 (£'), i.e. <S 1 (-E) is not invariant 
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Figure 2: (a) McGehee representation of the dynamics and the geometry 
of the manifolds in an energy surface H = E with E > near a saddle- 
centre equilibirium point. The concentric spheres represent pi — q\ = c 
(outer sphere), p\ — q\ = (dividing surface S^ S (E), middle sphere), and 
pi — qi = —c (inner sphere, not visible). The equator of the dividing surface 
is the periodic orbit S^ mM (E). W s f (E), Wf(E), W£{E), and W£(E) are the 
forward and backward branches of its stable and unstable cylinders (W£(E), 
and W?(E) are not visible), (b) Section of (a) with a plane of constant angle 
about the symmetry axis in (a). The arrows indicate the Hamiltonian vector 
field. Note that the vector field also has components out this the plane. (The 
energy surface H = E contains no equilibrium points.) Trajectories in the 
white regions correspond to non-reactive trajectories, (c) A deformation 
S% S (E) of the dividing surface S^ S (E) and its intersections with the stable 
and unstable manifolds of the Lyapunov periodic orbit S^ mM (E). The 
Hamiltonian vector field is tangent tg S% S (E) along a circle ^(E) marked 
by a dotted line (not completely visible), (d) The same section as in (b) 
now also showing S^ S (E) in this plane. The key point to note here is that 
^(E) (the intersection of the circle where the vector field is tangent to 
Sf s {E)) moves into the white, non-reactive region, (e) 0(E) divides S^ S (E) 
into two hemispheres D% s t (E) and D^ s h (E). (The extra "bubbles" shown 
in (e), that are not shown in (c), correspond to the intersection of 0(E) 



under the Hamiltonian vector field (S (E) is not a periodic orbit). It can be 
shown that if the deformation of S 2 S (E) is small enough then all points of 
S (E) have piqi < 0, i.e. they are contained in the complement of the two 
volumes enclosed by the forward and backward reactive spherical cylinders 
Wf(E) and Wb(E). Again intuitively this can be deduced from Figures Hfc)- 
(e), where the essential point is seen in (d) where S^ S (E) protrudes into the 
white, non-reactive region. This is where the trajectories become tangent 
to St s (E), i.e., at some point on S 2 S (E) trajectories "bounce off". 

Consider the part of S^ S (E) which is contained in the forward reactive 
spherical cylinder Wf(E) and marked by the light grey region in figure[2fc). 
The boundary of the light grey region corresponds to the intersection with 
the forward stable and unstable cylinder branches Wj(E) and WJ (E). The 
boundary is in general not a smooth manifold but it is homeomorphic to a 
circle S . It is possible to deform this piecewise smooth circle continuously 
onto the Lyapunov periodic orbit S^ mM (E) without leaving the stable and 
unstable cylinders. Stokes' theorem implies that the difference of the in- 
tegrals of <p along the piecewise smooth circle and the Lyapunov periodic 
orbit is given by the integral of £1' over the region on the stable and unstable 
manifolds swept out in the deformation process. The flux form Q' vanishes 
on the stable and unstable manifolds as its definition involves the interior 
product with the Hamiltonian vector field ( "there is no flux through invari- 
ant manifolds" ) . It thus follows that the flux through the light grey part of 
D 2 S { (E) in figure[2fd) is equal to the flux through D 2 g { (E). Therefore since 
the flux form f2' is strictly positive on the complete hemisphere D? f(E), the 
forward flux through S^ s (E) is larger than the forward flux through S 2 B (E) , 
and this completes the argument. 

All the arguments above immediately carry over to systems with more 
than 2 DOF by simply adjusting the dimension of the involved manifolds. 
The essential geometric conditions and relations amongst the manifolds all 
hold. 

The dividing surface of minimal forward flux is not unique, and this is 
the reason that we chose iso-energetic deformations that did not preserve 
the NHIM. For example the equation p\ = aqi with a (slightly) deviating 
from 1 again defines on the energy surface a (2n — 2)-sphere which is a 
dividing surface and which has the same forward flux as S^~ 2 {E). In our 
consideration we considered a generic deformation which is a deformation 
that changes the equator of S 2 ™~ 2 (E) such that it no longer coincides with 
the NHIM SnhimC^)- Non " 

generic deformation in this sense do not change 
the flux. Similarly, our transition state for 2 DOF coincides with a PODS 
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in general only along its equator, i.e. at the Lyapunov periodic orbit, and 
this is what matters for the flux. 



4 Algorithm for computing the flux 

Provided a generic non-resonance condition between the linear frequencies 
Uk is fulfilled the normal form Hamiltonian © assumes the simple form 
H(l, J2, . . . , J n ) where = {p\ + q%)/2, k = 2, . . . , n, are action variables 
associated with the bath coordinates. Like the saddle integral 2 the ac- 
tions Jfe are constants of the motion. Writing the flux form Q' in terms of 
action-angle variables we obtain the result that the forward flux through the 
transition state is given by 

T E = (27r) n - 1 V(S) 

where V(E) is the volume in the space of the bath actions (J2, . . . , J n ) en- 
closed by the contour H(0, J2, . . . , J n ) = E. In fact the flux can be inter- 
preted as the volume enclosed by a contour of constant energy E in phase 
space of a reduced system with one dimension less than the complete system. 
In terms of the normal form coordinates the reduced system is explicitly de- 
scribed by the invariant subsystem which has Qx = pi = 0. The normally 
hyperbolic invariant manifolds are just the energy surfaces of this reduced 
system which is referred to as activated complex in the chemistry literature. 
The dimensionless quantity .Fe/Zi™ -1 , where h is Planck's constant, is Weyl's 
approximation of the integrated density of states of the reduced system and 
can be interpreted as the number of "transition channels". 

In the linear case we have H = XX + X^fc=2 u kJk and the energy surface 
H = E encloses a simplex in (J2, . . . , J n ) whose volume leads to the well- 
known result (see e.g. |17j for an historical background) 

showing that the flux scales with E n ~ l for energies close to the saddle energy. 
The normal form allows to include the non-linear corrections to this result 
to any desired order. 

5 Conclusions 

In this letter we have given a geometrical argument showing that the flux 
across the dividing surface or transition state for n DOF systems constructed 
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in ^ [2] is a minimum, in the sense that the flux corresponding to forward 
reactive trajectories (which is equal in magnitude but opposite in sign to the 
flux corresponding to backward reactive trajectories) is a minimum. Hence 
it is the optimal dividing surface sought for by variational transition theory. 

The key point for the construction of the dividing surface is the existence 
of a normally hyperbolic invariant manifold (NHIM) which exists on the 
energy surface near a saddle-center-... -center equilibrium point and which 
mainly controls the dynamics nearby. The NHIM has the topology of a 
(2n — 3)-sphere which can be considered as the equator of the dividing 
surface which is a (2n — 2)-sphere. The NHIM divides the dividing surface 
into two "hemispheres" which are open (2n — 2)-discs. The hemispheres are 
transverse to the Hamiltonian flow. Hence the dividing surface is everywhere 
transverse to the Hamiltonian flow except for its equator (the NHIM) which 
is an invariant manifold. The NHIM can be considered as the energy surface 
of an invariant subsystem ( "activated complex" ) with one DOF less than the 
complete system. The flux is the phase space volume enclosed by the energy 
surface of this reduced system. 

The dividing surface is not unique. Any (2n — 2)-sphere which con- 
tains the NHIM as its equator and which is transverse to the Hamiltonian 
flow except for its equator qualifies for a dividing surface. Nevertheless, 
all these dividing surfaces lead to the same flux. This makes it difficult (if 
not impossible) to compute a dividing surface from a variational principle 
of codimension 1 manifolds of the energy surface. Similarly, the NHIM ap- 
pears only as an indefinite critical "point" of the variational principle of 
codimension 2 manifolds of the energy surface of MacKay |18j what makes 
it practically infeasible to compute the NHIM from a variational principle. 
The normal form approach of ^ [2] is currently the only method to deter- 
mine these manifolds and also to compute the flux for which we gave an 
algorithm in this letter. 
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